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Using perturbative Quantum Chromodynamics, we compute dihadron fragmentation functions 
for a large invariant mass of the dihadron pair. The main focus is on the interference fragmentation 
function Hf, which plays an important role in spin physics of the nucleon. Our calculation also 
reveals that Hf and the Collins fragmentation function have a closely related underlying dynamics. 
By considering semi-inclusive deep-inelastic scattering, we further show that coUinear factorization 
in terms of dihadron fragmentation functions, and collinear factorization in terms of single hadron 
fragmentation functions provide the same result in the region of intermediate invariant mass. 



PACS numbers: 12.38.Bx, 12. 39. St, 13.85.Ni, 13.87.Fh 

/. Introduction. — Fragmentation functions (FFs) for 
quarks and gluons parameterize the hadronization tak- 
ing place in high-energy scattering processes with identi- 
fied strongly interacting particles in the final state. The 
main focus is typically on FFs describing the transition 
of a parton into a single hadron — see Ref. [l[ for a field- 
theoretic definition of such objects. However, already in 
the late 1970's dihadron fragmentation functions (DiFFs) 
were introduced in order to quantify the hadron structure 
of jets Moreover, it was shown that DiFFs are needed 
to obtain a consistent result for the production of two 
hadrons in electron-positron annihilation when working 
beyond leading order in perturbative Quantum Chromo- 
dynamics (QCD) 3]. In the meantime, DiFFs also play 
a considerable role in heavy ion physics — see Q and 
references therein. 

In 1993, it was proposed 0] that quark fragmentation 
into two hadrons can also be used to address the transver- 
sity distribution hi of the nucleon 0, Q- To this end, 
one can study the production of two hadrons in semi- 
inclusive deep-inelastic scattering (DIS) in the current 
fragmentation region . If the target is transversely po- 
larized, there exists a correlation between the spin vector 
of the target and the orientation of the plane given by 
the momenta of the two hadrons. This observable con- 
tains the product of hi and a new fragmentation func- 
tion {Hf in the notation of Ref. [1]), which describes the 
strength of a correlation between the transverse polariza- 
tion of the fragmenting quark and the orientation of the 
hadron plane. Like the transversity, is chiral-odd, 
and it results from the interference between two different 
production amplitudes why it is normally referred to as 
interference fragmentation function in the literature j9| . 

Data on the mentioned observable in semi-inclusive 
DIS have already been obtained by the HERMES and 
COMPASS Collaborations \^^. The major difficulty 
is that, a priori^oth. hi and are unknown. Existing 



models for 

are still in a too early stage for 
getting a quantitative constraint on the transversity. (For 
a related discussion we refer to However, as was 

shown in [l5j , one can measure two back-to-back hadron 



pairs in electron-positron annihilation in order to get a 
handle on H? — see also [H, [13, and Ref. (3 for pre- 
liminary data from the Belle Collaboration. A combined 
analysis of the two processes allows one, in principle, to 
extract both unknown functions. (First steps towards 
this goal are outlined in [l^.) Such a strategy would 
be very similar to the combined analysis of the Collins 
effect [20| in semi-inclusive DIS and in electron-positron 
annihilation 



21 



from which first information about 



the transversity was obtained |22j . 

The interference FF Hi and the Collins function [Hi 
in the notation of Refs. 23, [l^l) can be considered as 



complementary tools for getting a handle on the transver- 
sity distribution, with both having advantages and draw- 
backs. An important advantage in the case of is 
the fact that one can integrate over the total trans- 
verse momentum of the two hadrons in the final state, 
leading to a collinear factorization formula. In con- 
trast, the Collins effect relies on factorization in terms 
of transverse momentum dependent parton correlators 
(TMD-factorization) [isl - fiit . which has additional tech- 
nical complications. On the other hand, when using Hi , 
the dependence on the relative transverse momentum of 
the two hadrons must be kept in order to have a well- 
defined hadron plane. This implies that Hf must also 
depend on the invariant mass M^i^ of the dihadron sys- 
tem. If M}ih is of the order of h.QCD , DiFFs are entirely 
non-perturbative objects. In this kinematical region, one 
can only fit the DiFFs to experimental data or try to 
estimate them by using some model for the strong inter- 
action in the non-perturbative regime [l^, [l^l ■ (As a 
matter of principle, FFs cannot be computed in lattice 
gauge theory.) 

In the present paper, we apply perturbative QCD in or- 
der to evaluate DiFFs for Af^^ :s> Aqcd- The main focus 
is on the interference FF H^. For large Mf^h, the DiFFs 
can be expressed as a convolution of hard coefficients 
and (collinear) single hadron fragmentation correlators. 
In particular, the calculation determines the behavior of 
DiFFs as a function of Mhh- While the unpolarized DiFF 
Di drops like the interference DiFF H^ behaves 
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FIG. 1: Graphical representation of dihadron fragmentation 
function for a quark. 



like 1 /M^^ . We also argue that and the Collins func- 
tion Hi (at large transverse momentum 29]) depend on 
the same two coUinear twist-3 fragmentation correlators, 
showing that the underlying dynamics of both functions 
is closely related. In addition, we compute the (trans- 
verse) spin dependent cross section for dihadron produc- 
tion in semi-inclusive DIS for Aqcd ^ ^ Q {Q de- 
noting the virtuality of the exchanged photon) with single 
hadron FFs. By comparing this result with the cross sec- 
tion obtained in the framework of DiFFs, we show that 
coUinear factorization in terms of DiFFs holds as long as 
Mhh <C Q. In its spirit, our study is similar to recent 
work in which certain transverse momentum dependent 
parton correlators were evaluated for large transverse mo- 
menta, and the matching between coUinear factorization 
and TMD-factorization was explicitly shown for interme- 



see also 



for an 



diate transverse momenta 
overview. 

//. Kinematics and definition of dihadron fragmentation 
functions. — We start by discussing the kinematics for 
the fragmentation of a quark into two hadrons (displayed 



for the squared amplitude also in Fig.[T]), 

qik) ^ hi{Pih) + h2{P2h) + X . 



(1) 



We assume that the quark has a large light-cone minus 
momentum . For later convenience, we choose a ref- 
erence frame in which one hadron has no transverse mo- 
mentum. In the hadron-1 frame, for instance, the light- 
cone components of the hadron momenta can be repre- 
sented as 



Pih = 

P2h = 



0, zihk , 



ML 
2zihk- 



Z2hk , X —Mhh 

V Zlh 



(2) 



Neglecting the hadron masses, one readily verifies that 
{Pih + P2h)'^ = We also introduce the total hadron 

momentum as well as the momentum difference according 
to 



P} 



hh 
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Ih 



P 



2h ■ 



R 



Pih 



2h 



(3) 



Their minus momenta are given by P^^^ = Zhhk , R 
Zhhk~, with Zhh = zih + Z2h and Zhh = [zih - Z2h)/2. 



The operator definition of the unpolarized DiFF and 
the interference DiFF, for a quark fiavor q, reads 



hh 



^hh 



^Ih X 



2tt 



32(27r)3 

X {0\i;'^{y+)\Pih,P2h,X){Pih,P2h,X\tlj''{0)\Q) 



l^D\{zhh,ML) 



Rra 



2\Rt\ 



Ht\zhh,ML), (4) 



where a gauge link has been suppressed. Note that our 
definition of differs from the one in Ref. [8] by some 
prefactors. We have integrated over the relative longi- 
tudinal momentum fraction Zhhi which is needed if one 
wants to apply the coUinear factorization discussed in the 
next section. 

With these conventions, the parton model cross section 
for the production of two hadrons in semi-inclusive DIS 
(with a transversely polarized proton), ep^ — >■ ehih2X, 
takes the form 



da 



dxB dyd(j)s dzhhdM^ d<j)R 



2al^sxB 2 
r)4 Z^^i 



1-2/ 



f1{xB)Dl{zhh,MD 



(1 - 2/)sin(0fl + 4>s)h\{xB)Ht''{zhh,MD 



(5) 



Equation ([S]) contains both the unpolarized cross section 
and one component depending on the transverse target 
polarization. The azimuthal angle between Rt and the 
lepton plane is denoted by (pR , the azimuthal angle of the 
transverse spin vector of the proton is denoted by 0s, 
while XB and y are the commonly used DIS variables. 
///. Dihadron fragmentation functions at large invariant 
mass. — When Mhh ^ ^qcd, dihadron fragmentation 
of a quark can be viewed as a two-step process: first, the 
quark splits into a quark (with momentum Iq) and a gluon 
(with momentum Ig), which is calculable in perturbative 
QCD. Second, each of these two partons fragments into 
a single hadron. This scenario is illustrated in Fig. [5] 
We introduce momentum fractions (zi,Z2) through 



In 



Pih 
Z\ 



In 



Pih 
Z2 



(6) 



and define ^ = zih/zi. Evaluating the diagram in 
Fig.HKa) one finds for the unpolarized DiFF 



Dl{zhh,MD = 



dz 



d^ 



(7) 



I D^^^^ representing unpol; 
hadron FFs. A second term, where hadron 1 originates 



with d'i^^'' and D^^^^ representing unpolarized single- 
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FIG. 2: Generic diagrams contributing to DiFFs for large in- 
variant mass. For the unpolarized DiFF Di the diagram (a) 
needs to be considered, whereas in the case of Hi both dia- 
grams are relevant. 



from the fragmentation of the gluon, is not written for 
brevity. The result in Eq. dH) shows that the DiFF Di 
drops like 1 /M^j^ . Note also that (O allows one to recover 
the inhomogeneous part of the evolution equation for Di 



by integrating over M^^ 

Computing for large Mhh is much more involved as 
one has to use collinear twist-3 factorization [3^[33|. One 
contribution arises from the diagram in Fig.l^Ja), if the 
relative transverse momentum between the quark with 
momentum Ig and the hadron 1 is kept. Also, 3-parton 
correlators including a transverse gluon field A± enter 
the calculation — see the sample diagram in Fig.^b). 
The correlators associated with these two contributions 
(the so-called 9^-contribution and the j4^-contribution) 
are of the form {'ipd^'ip) and ■(/)), respectively. In 

this paper, contributions from 3-gluon correlators are ne- 
glected since they do not affect any of our general results. 
We leave this part, details of the calculation, and the dis- 
cussion about a potential singularity in ([T]) at ^ = and 
^ = 1 for future work 35 1. 

The calculation of the hard coefhcient for the A^^- 
contribution is essentially identical to the corresponding 
part in the treatment of the Collins function at large 
transverse momentum [29j. On the other hand, to ob- 
tain the c?j_-contribution is more complicated than for 
the Collins function evaluation. To this end, we assign a 
(relative) transverse quark momentum according to 



P' 



9-1- ' 



(8) 



where z'^ — zi+Szi, and P{i^ = Pih + SPih- To keep Mhh 
fixed, P2h must also change, i.e., = + SP: 



2h- 



In 

particular, 5P2h± = —Z2lq±. The kinematics is entirely 
determined by the constraints 



Pih ■ SPih 

Plh ■ SP2h 



0, 

P2h 



SP2k = 



P2h ■ SP2h 
SPlh - , 

SPih± = , 



0, 



(9) 



where we use the on-shell condition for the two hadrons, 
constraints from keeping Mhh and Zhh fixed, and the fact 
that we are working in the hadron- 1 frame. (Recently, 
we used a related approach for computing a particular 



single spin asymmetry in the Drell-Yan process 36[.) The 
solution to the set of equations in ([9)) reads 



6P 



2h 



P2h± Zhh 
P2h±lq± Z2 _'2k^lq±Z2ZlhZ2h 
Zhh P2h± Zhh 



5ZI = — Z1Z2(Z2 - Zl)(l - C)4^ 
Zhh r2hA 



, -Z2lq^j , 

(10) 



The partonic scattering amplitude M depends now on 
lq±. We expand M, 



M{P[^^,Pf,^,z[) = M{P,h,P2h,z,) 



dMiPih,Pk,z[) 



dl„ 



lq± + 



(11) 



la±=0 



M 35 



for obtaining the relevant 



and keep the term linear iiiJ 
(twist-3) 9j_-contribution 

Both the 9^-contribution and the ^^-contribution can 
be brought into a gauge invariant form. After collecting 
all the pieces we find [35| 



H^''izhh,M^f,) 



dzhh 



2ttM^^ J_1 ^ZlhZ2h 7^1, C 



1-Z2h 



xA'^^/^Zhh,OD'i'^'{^), (12) 
where the function A is defined as 



A = Cf 



3 d Hizi)\ ^zi ~ Z2 



dzi zf 



dzi 



PV 



Zhh 
Hf{z1:Zi) 



H{zi) 



2e_ 

1-C 



Zl V Zl Zl 

Ca 2zih zizijzi + zi) - zihjzf + zf) 
2 Zl zi{zi - zi){zi - zih) 



(13) 



This result shows that behaves like 1/M|^ for large 
Mhh- It also reveals an intimate connection between 
and the Collins function (at large transverse momen- 
tum) [2^ , as both functions depend on the same collinear 
twist-3 correlation functions H and Hp, which we take 
in the definition of Ref. [2^ . The hard coefficients (of the 
9_L-contribution) differ in both cases. 
IV. Single spin asymmetry for dihadron production in 
semi-inclusive DIS — In this section, we investigate the 
validity of the factorization formula (5) for the cross sec- 
tion of dihadron production in semi-inclusive DIS for the 
case Aqcd ^ Mhh ^ Q- Though one might expect 
this factorization to hold, to the best of our knowledge 
no explicit supportive calculation exists. We focus on 
the discussion of the spin-dependent component of the 
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FIG. 3; Sample diagrams contributing to dihadron produc- 
tion in semi-inclusive DIS. Diagram (a) generates a leading 
power contribution and diagram (b), in the light-cone gauge, 
is power-suppressed for Mhh <^ Q- 



cross section, which is nontrivial already at lowest order 
in perturbation theory. 

For Mhh of the order of Q, one can use collincar fac- 
torization in terms of single hadron fragmentation corre- 
lators. Sample diagrams are shown in Fig. 3. We have 
evaluated the cross section for this kinematics, and then 
expanded the result for Aqcd ^ Mhh ^ Q- For factor- 
ization in terms of DiFFs to hold, the result has to match 
with (5), if iJj^ for Mhh > ^qcd from ((T2|) is inserted. 
By inspecting the Feynman diagrams it becomes obvious 
that this matching is indeed nontrivial. For instance, di- 
agram (b) in Fig. 3 has no counterpart in the calculation 
of for large Mhh- However, it turns out that for 
Mhh ^ Q this diagram is power-suppressed. This is true 
only in the light-cone gauge — which we use for 
this analysis. In a covariant gauge, the treatment gets 
more involved [35| . In the end we indeed find a match- 
ing, showing the consistency of the factorization in terms 
of DiFFs for Mhh > Aqcd as long as Mhh < Q- For 
the unpolarized cross section the corresponding analysis 
is trivial to lowest order, but also becomes nontrivial ones 
loop corrections are included. 

IV. Conclusions. — In this paper, we have studied 
DiFFs for a large invariant mass Mhh of the dihadron 
pair, where perturbative QCD can be applied. The 
main focus has been on the interference DiFF , which 
drops like 1/M^f^ for large Mhh and is related to the 
same universal twist-3 collinear fragmentation correla- 
tors that describe the Collins FF (at large transverse 
momentum). The analysis also predicts that the trans- 
verse single spin asymmetry for dihadron production in 
semi-inclusive DIS behaves like 1/Mhh- The preliminary 



COMPASS data [ll|, ranging up to Mhh « 2GeV, are 



in agreement with this general result. We expect a corre- 
sponding behavior for the so-called Artru-Collins asym- 
metry in electron-positron annihilation 15|, for which 
preliminary data from Belle exist Is*]. For the case of 
semi-inclusive DIS we have also shown explicitly, to low- 
est nontrivial order in perturbation theory, that collinear 
factorization in terms of DiFFs is consistent for large Mhh 
provided that Mhh Q- 
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